Search for: All records

Abstract Observations of quasar absorption spectra provide strong evidence that reionization extended belowz= 6. The relationship between Lyαforest opacity and local galaxy density (the opacity-density relation) is a key observational test of this scenario. Using narrow-band surveys ofz≈ 5.7 Lyαemitters (LAEs) centered on quasar sight lines, ref. [1] showed that two of the most transmissive Lyαforest segments at this redshift intersect under-densities in the galaxy distribution. This result is in tension with models of a strongly fluctuating ionizing background, including some models of late reionization, which predict that the vast majority of these segments should intersect over-densities where the ionizing intensity is strongest. In this paper, we use radiative transfer simulations to explore in more detail the opacity-density relation predicted by late reionization models. We find that fields like the one toward quasar PSO J359-06 — the more under-dense of the two transmissive sight lines in ref. [1] — are typically associated with recently reionized gas inside of cosmic voids where the hotter temperatures and rarefied densities enhance Lyαtransmission. The opacity-density relation's transmissive end is sensitive to the amount of neutral gas in the voids, as well as its morphology, set by the clustering of reionization sources. These effects are, however, largely degenerate with each other. We demonstrate that models with very different source clustering can nonetheless yield nearly identical opacity-density relations when their reionization histories are calibrated to match Lyα forest mean flux measurements atz< 6. In models with fixed source clustering, a lower neutral fraction increases the likelihood of intersecting hot, recently reionized gas in the voids, increasing the likelihood of observing fields like PSO J359-06. For instance, the probability of observing this field is 15% in a model with neutral fractionx_HI= 5% atz= 5.7, three times more likely than in a model withx_HI= 15%. The opacity-density relation may thus provide a complementary probe of reionization's tail end. 

Abstract The wealth of high-quality observational data from the epoch of reionization that will become available in the next decade motivates further development of modeling techniques for their interpretation. Among the key challenges in modeling reionization are (1) its multi-scale nature, (2) the computational demands of solving the radiative transfer (RT) equation, and (3) the large size of reionization's parameter space. In this paper, we present and validate a new RT code designed to confront these challenges.FlexRT(Flexible Radiative Transfer) combines adaptive ray tracing with a highly flexible treatment of the intergalactic ionizing opacity. This gives the user control over how the intergalactic medium (IGM) is modeled, and provides a way to reduce the computational cost of aFlexRTsimulation by orders of magnitude while still accounting for small-scale IGM physics. Alternatively, the user may increase the angular and spatial resolution of the algorithm to run a more traditional reionization simulation.FlexRThas already been used in several contexts, including simulations of the Lyman-αforest of high-zquasars, the redshifted 21cm signal from reionization, as well as in higher resolution reionization simulations in smaller volumes. In this work, we motivate and describe the code, and validate it against a set of standard test problems from the Cosmological Radiative Transfer Comparison Project. We find thatFlexRTis in broad agreement with a number of existing RT codes in all of these tests. Lastly, we compareFlexRTto an existing adaptive ray tracing code to validateFlexRTin a cosmological reionization simulation. 

Abstract During reionization, intergalactic ionization fronts (I-fronts) are sources of Lyαline radiation produced by collisional excitation of hydrogen atoms within the fronts. In principle, detecting this emission could provide direct evidence for a reionizing intergalactic medium (IGM). In this paper, we use a suite of high-resolution one-dimensional radiative transfer simulations run on cosmological density fields to quantify the parameter space of I-front Lyαemission. We find that the Lyαproduction efficiency — the ratio of emitted Lyαflux to incident ionizing flux driving the front — depends mainly on the I-front speed and the spectral index of the ionizing radiation. IGM density fluctuations on scales smaller than the typical I-front width produce scatter in the efficiency, but they do not significantly boost its mean value. The Lyαflux emitted by an I-front is largest if 3 conditions are met simultaneously: (1) the incident ionizing flux is large; (2) the incident spectrum is hard, consisting of more energetic photons; (3) the I-front is traveling through a cosmological over-density, which causes it to propagate more slowly. We present a convenient parameterization of the efficiency in terms of I-front speed and incident spectral index. We make these results publicly available as an interpolation table and we provide a simple fitting function for a representative ionizing background spectrum. Our results can be applied as a sub-grid model for I-front Lyα emissions in reionization simulations with spatial and/or temporal resolutions too coarse to resolve I-front structure. In a companion paper, we use our results to explore the possibility of directly imaging Lyαemission around neutral islands during the last phases of reionization. 

Abstract Long troughs observed in thez> 5.5 Lyαand Lyβforests are thought to be caused by the last remaining neutral patches during the end phases of reionization — termed neutral islands. If this is true, then the longest troughs mark locations where we are most likely to observe the reionizing intergalactic medium (IGM). A key feature of the neutral islands is that they are bounded by ionization fronts (I-fronts) which emit Lyman series lines. In this paper, we explore the possibility of directly imaging the outline of neutral islands with a narrowband survey targeting Lyα. In a companion paper, we quantified the intensity of I-front Lyαemissions during reionization and its dependence on the spectrum of incident ionizing radiation and I-front speed. Here we apply those results to reionization simulations to model the emissions from neutral islands. We find that neutral islands would appear as diffuse structures that are tens of comoving Mpc across, with surface brightnesses in the range ≈ 1 - 5× 10^-21erg s^-1cm^-2arcsec^-2. The islands are brighter if the spectrum of ionizing radiation driving the I-fronts is harder, and/or if the I-fronts are moving faster. We develop mock observations for current and futuristic observatories and find that, while extremely challenging, detecting neutral islands is potentially within reach of an ambitious observing program with wide-field narrowband imaging. Our results demonstrate the potentially high impact of low-surface brightness observations for studying reionization. 

Abstract Recently, the mean free path of ionizing photons in the z = 6 intergalactic medium (IGM) was measured to be very short, presenting a challenge to existing reionization models. At face value, the measurement can be interpreted as evidence that the IGM clumps on scales M ≲ 10 8 M ⊙ , a key but largely untested prediction of the cold dark matter (CDM) paradigm. Motivated by this possibility, we study the role that the underlying dark matter cosmology plays in setting the z > 5 mean free path. We use two classes of models to contrast against the standard CDM prediction: (1) thermal relic warm dark matter (WDM), representing models with suppressed small-scale power; (2) an ultralight axion exhibiting a white noise-like power enhancement. Differences in the mean free path between the WDM and CDM models are subdued by pressure smoothing and the possible contribution of neutral islands to the IGM opacity. For example, comparing late reionization scenarios with a fixed volume-weighted mean neutral fraction of 20% at z = 6, the mean free path is 19 (45)% longer in a WDM model with m x = 3 (1) keV. The enhanced power in the axion-like model produces better agreement with the short mean free path measured at z = 6. However, drawing robust conclusions about cosmology is hampered by large uncertainties in the reionization process, extragalactic ionizing background, and thermal history of the Universe. This work highlights some key open questions about the IGM opacity during reionization.